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ABSTRACT 

Infrared absorption lines of }i^ , including the metastable i?(3,3)' line, have been observed toward 
eight bright infrared sources associated with hot and massive stars located in and between the Galactic 
Center Cluster and the Quintuplet Cluster 30 pc to the east. The absorption lines with high velocity 
dispersion arise in the Galaxy's Central Molecular Zone (CMZ) as well as in foreground spiral arms. 
The temperature and density of the gas in the CMZ, as determined from the relative strengths of 
the H;^ lines, are T=200 - 300 K and n <50 - 200 cm~'^. The detection of high column densities 
of Hg toward all eight stars implies that this warm and diffuse gaseous environment is widespread 
in the CMZ. The products of the ionization rate and path length for these sight lines are 1000 and 
10 times higher than in dense and diffuse clouds in the Galactic disk, respectively, indicating that 
the ionization rate, C, is not less than 10~^^ s~^ and that L is at least on the order of 50 pc. The 
warm and diffuse gas is an important component of the CMZ, in addition to the three previously 
known gaseous environments: (1) cold molecular clouds observed by radio emission of CO and other 
molecules, (2) hot (T = 10^ - 10^ K) and highly ionized diffuse gas (rig = 10 - 100 cm~'^) seen 
in radio recombination lines, far infrared atomic lines, and radio-wave scattering, and (3) ultra-hot 
(T = 10^ - 10^ K) X-ray emitting plasma. Its prevalence significantly changes the understanding of 
the environment of the CMZ. The sight line toward GC IRS 3 is unique in showing an additional Hj^ 
absorption component, which is interpreted as due to either a cloud associated with circumnuclear disk 
or the "50 km s^^ cloud" known from radio observations. An infrared pumping scheme is examined 
as a mechanism to populate the (3,3) metastable level in this cloud. 
Subject headings: Galaxy: center — ISM: clouds — ISM: lines and bands — ISM: molecules — stars: 

individual (GC IRS IW, GC IRS 21, GC IRS 3, NHS 21, NHS 22, NHS 42, NHS 25, 

GCS 3-2) 



1. INTRODUCTION 

1.1. The Central Molecular Zone 

The region of the Galaxy within 200 pc of the cen- 
ter, known as the Central Molecular Zone (CMZ), con- 
tains 10% of the Galactic interstellar molecular mass 
(iMorris fc Serabvn! il996: Gcnzcl. HoUenbach fc Townes I 
11994 iMezger Duschl. fc Zvlka nJOel . and references 
therein). The interstellar medium (ISM) in the CMZ 
is exceptional in many other ways. It is highly turbu- 
lent, pervaded by an intense magnetic field (for a re- 
view; 'Morris 2006) , contains a multitu de of massiv e stars 
e.g. fMoneti et al..il992l:ICotera et al.lll998.: .Muno et al 



2006) 



199C 



and super-massi ve clusters (e.g.. iNagata et ah 
, 119951: iFigeil [200l . for a short review), and also 
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is home to numerous highly energetic objects rang- 
ing from s upernova remnants to cataclysmic variables 
(CVs, e.g., iMuno et aDl2003f ). The particle densities 
in its dense molecular clouds are of order lO'* cm~'^ 
and higher. The volume filling factor, /v, of dense 
mole cular clouds has previous ly been estimated to be ~ 
10% (iMorris fc Serabvnl [l996D. but is more likely ~ 1% 
(|Oka et al.ll2005f) in view o f the observed visu al extinc- 
tion of ^y = 25 - 40 mag (jCotera et al.l |200(t). 

If dense molecular clouds only occupy 1% of the CMZ, 
what interstellar environments make up the remaining 
99%? There have been some reports of lower density 
"diffuse" molecular gas in the intercloud region. Detailed 
analysis of the J = 1 — > e mission lines of CO and 
Ci^O bv lDahmeneTall (fT99g). and the J = 2 -> 1 and 
J = 1 ^ emission hues of CO by Oka et all (|1998af) 
have demonstrated the presence of gas w ith densities of 
~ 10^-^ cm-3. Indeed lOka et all U2M^ conclude that 
CO emission from the Galactic center (GC) is dominated 
by emission from gas of that densit y. The volume fillin; 
factor of /v ^ 0.1 claimed bv Morris fc Serabvnl (jl99i 
may include such gas. In this paper we present evidence 
for an even lower density molecular environment in which 
H2 must still be plentiful but where CO is not abundant. 
In hindsight, the radio absorpt ion lines with h igh velocity 
dispe rsion of OH (e g. .Bolton et all I1964D. H9CO (e 
g. P almer et al.lll969D, HCO+ (iLinke. Stark, fc FrerkingI 
119811), and CO (|Oka et al.i ri998b|) mav originate in this 
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newly reported environment. 

A second candidate to fill the intercloud medium is the 
ultra-hot plasma (T -- lO"^ - 10^ K) observed by its dif- 
fuse X-ray emission at energies of 0.5 - 10 keV. The radi- 
ation from the plasma is observed in all sightlincs to the 
CMZ (e.g.. Kovama et al. 1989: Yamauchi fc Kovamal 
119931 iKovama et al.lll996[ r iLazio fc CordesI (|1998f) sug- 
gested that the highly ionized regions with Ug ~ 10 cm~'^ 
and T r^ 10® K that are required to explain the A^ depen- 
dent radio- wave scattering toward the GC are located at 
the photon-dominated interfaces between the molecular 
clouds and this ultra-hot X-ray emitting plasma. They 
proposed that high density molecular clouds with a vol- 
ume filling factor of perhaps /v ^ 0.1 are surrounded by 
the scattering electron gas and the rest of the space is 
filled with this plasma (see their Fig. 9). 

For an environment such as that described above, 
which is so dominated by ionized gas, "Central Plasma 
Zone" rather than CMZ would be a more appropriate 
term for the region. For a variety of reasons the ultra-hot 
and hot plasma cannot occupy the same regions where 
molecules abound. How those different categories of gas 
are distributed and how they coexist in the CMZ is yet 
to be understood. Observations of a new astrophysical 
probe of the CMZ, the infrared absorption spectrum of 
Hj]" ([Geballc & Oka 199 a), are beginning to shed fresh 
hght on these questions (jBoldvrev fc Yusef-Za deh 2006) . 
Earlier measu rements of this molec ul ar ion in the Galac- 
tic center by ^eballe et all (|1999( l. iGoto et al.l (|2002f ). 
and lOka et al.l ( 20051 ) were made along only a very few 
sight lines. In this paper we report and analyze observa- 
tions of seven transitions of H;^ along eight sight lines, 
which provide a more comprehensive diagnostic of the ex- 
tent and nature of the H^ - containing gas in the CMZ. 

1.2. H^ as an Astrophysical Probe 

The crucial information obtained from the observations 
of Hjj' reported here is the product of the ionization rate 
of molecular hydrogen (, and the aggregate length of the 
line of sight through the Hj^ - containing clouds, L. The 
simple chemistry of Hg allows us to obtain this product 
from the observed total Hj^ column densities (lOka et al.l 
120051 IOkall2006bl ). as described below, under the as- 
sumptions that (1) there is a steady-state balance be- 
tween creation and destruction of H^^ and (2) the domi- 
nant supplier of electrons is carbon, all of which is singly 
ionized. 

Creation of H^ begins with ionization of H2. Be- 
cause in interstellar regions where H2 is plentiful the ion- 
neutral reaction H2 + H^ -^ H^ + H is rapid, the limit 
to the formation of Hg is the much slower ionization rate 
of molecular hydrogen, 

H2-fR 



Hj 



R'. 



The ionizing rays, R, are dominated by cosmic rays in 
the Galactic disk ()Indriolo et al.l l2007l and references 
therein), but in the CMZ the ionization by X-rays and 
EUV radiation may be significant. The production rate 
per unit volume of Hg" is thus given by ^n(H2), where 
C is the effective ionization rate of H2 including all the 
aforementioned contributors. In steady state, produc- 
tion of H;^ balances its destruction, which is dominated 
by dissociative recombination with electrons. 



H+-he^3H or H2 -H H, 
and one obtains (jGeballe et al.|[l999l ) 

Cn(H2) = fcon(H+)n(e), 

where ke is the electron recombination rate constant, 
whose temperature dependent value is ~ 10^^ cm'^ s~^ 
(jMcCall et al.ll20M ). Replacing the electron density with 
the atomic carbon density after depletion gives 

n(e)/n(H2) = (2//) [nc/nn]syRx, 

where Rx = 3 - 10 is the factor increase in relative 
carbon abundance at t he Galactic center ov er the so- 
lar vicinity, [nc/'^Hlsv (jSodroski et al.lfl995f ). and / = 
2n(H2)/nH is the fractional density of molecular hydro- 
gen relative to the total number of hydrogen atoms, 
riH = 2n(H2) -l-n(H). Using iV( H+) = £n(H +) because 
of the constancy of n(R^) (e.g.. IOka1l2006bl ). we obtain 
the equation central to the analysis in this paper relating 
the product (^L to the total H^ column density and other 
measurables, 

(L = 2fce hc/nnlsv RxN{R+)totai/f. (1) 

As discussed bv lOka et al.l ()2005[ ). in the Galactic cen- 
ter the spectrum of Hg provides direct information about 
the temperature and density of the gas. There the rela- 
tive population of H;^ in the (J, K) — (3,3) metastable 
level, which is 361 K above the lowest (1,1) level (note 
that J — does not exist in the ground vibrational 
state), is measurable and is a useful thermometer for 
temperatures of ^ 100 - 1000 K. Observations of the frac- 
tional population in the (2,2) unstable level, from which 
R^ decays to the (1,1) level by spontaneous emission 
in 27 days (|Pan fc Okalll986l : iNeale. Miller, fc TennvsoiTI 
I1996D . provides a good measure of cloud densities less 
than several hundred per cm^. The lifetime of Hg in 
diffuse clouds is l/{konc) « 10^ s, which is two orders 
of magnitude longer than the collisional timescale with 
molecular hydrogen. The rotational temperature of Hg 
is therefore maintained via collisions with the ambient 
H2. The thermalization of H^ in interstellar clouds has 
been modeled by lOka fc EppI ()2004l ). who provide H^ 
population ratios n(3,3)/n(l,l) and n(3,3)/n(2,2) plot- 
ted as functions of temperature and density (see their 
Figs. 2-4). 

2. OBSERVATIONS 

Bright infrared sources corresponding to dust shells 
surrounding luminous stars or to hot and bright stars 
with a dearth of photospheric abso rption features were 
selected from the Sgr A* cluster (" Beckhn et al.l 119781: 
iKrabbe et anfl99 5: Tanner et al. 20 0J), the Q uintuplet 
cluster (|Okuda et al.lll99 0: Naga ta et al.ll990D . and can- 
didates in their vicinity (Nagata et al.l Il993l ) to use as 
probes of foreground Hj^ (Table [1]). Spectroscopic ob- 
servations were carried out at the Subaru Telescope over 
several nights between 2004 July 7 UT and 2004 Sep 26 
UT, using the Infrared Camera and Spectrogr aph (IRCS) 
(jTokunaga et al.lll998l:lKobavashi et al.ll2000D . The data 
obtained on 2001 Jun 16 UT, originally published in 
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Source 



TABLE 1 
Target list. 



R.A. 

(J2000) 



Dec. 

(J2000) 



L/L' Other Names 
[mag] in SIMBAD 



Reference 
1: 2, 3. 



GC IRS 21 . . . 17:45:40.2 -29:00:31 -0.0561 -0.0468 NAME SGR A IRS 21, GCIRS 21 1, 2, 3. 

GCIRS3... 17:45:39.8 -29:00:24 -0.0552 -0.0448 5.3 BHA 11, GCIRS 3 1,2,3. 

GC IRS IW . . . 17:45:40.4 -29:00:27 -0.0549 -0.0471 5.5 NAME SGR A IRS IW, GCIRS IW 1, 2, 3. 

NHS 21... 17:42:53.8 -28:51:41 +0.0992 -0.0553 4.62 qF 577, [NHS93] 21 4,5,6. 

NHS 22... 17:46:05.6 -28:51:41 +0.1199 -0.0482 6.0 [NHS93] 22 4,5,6. 

NHS 42... 17:46:08.3 -28:49:55 +0.1481 -0.0426 6.05 qF 578, [NHS93] 42 4,5,6. 

NHS 25... 17:46:15.3 -28:50:04 +0.1592 -0.0657 6.2 NAME PISTOL STAR, [NHS93] 25 4,5,6. 

GCS 3-2 .. . 17:46:14.7 -28:49:41 +0.1636 -0.0606 2.5 GCS 3-2, [NHS93] 24 4, 5, 6. 



Not e. — Reference: 1. IKrabbe et al"] l|1995l) . 2. lOtt, Eckart, &: Genzell II1999D . 3. ITanner et~all ll2006f) . 4. INagata et al.l J1990I1 . 5. INagata et al.l 
119931) . 6. IFJger. McLean, fc Morrij 1119991) ! 



I Goto et al] (|2002f ). were combined with the newer data. 
A summary of the observations is given in Table [21 

A curvature sensing adaptive optics (AO) system 
(jTakami et al.i r2004) was used whenever a wavefront ref- 
erence star was available. Use of the AO system signifi- 
cantly increased the throughput of starlight through the 
narrow (C/IS) slit, which was necessary to attain high 
spectral resolution (A/AA = 20,000; Av = 15 km s"^). 
In order to subtract sky and telescope emission, either 
the telescope itself or the tip-tilt mirror in the AO sys- 
tem was nodded by about 3" along the slit so that the 
stellar spectrum fell on two different parts of the camera. 
The slit position was set along the east-west axis in all 
cases except when nearby faint sources fell onto the slit, 
which would have hampered a clean sky-subtraction. 

The angle of cross-dispersing grating in IRCS was ad- 
justed so that as many Hj lines as possible were observed 
in the time available. The major absorption lines essen- 
tial for the gas diagnostics (the i?(l,0), i?(l,l)", i?(l,l)', 
(3(1,0), and (3(1,1) lines from the lowest (1,1) and (1,0) 
levels, and the i?(3,3)' and i?(2,2)' lines from the higher 
(3,3) and (2,2) levels) were covered in two grating set- 
tings. The (5(1,0) transition at 3.9530 fiui was covered 
by both of the gr ating settings. De tails of the spectral 
setup are found in lGoto et al.l (|2002D . 

Standard stars with early spectral types were observed 
on the same nights through similar airmasses as the 
Galactic center stars in order to remove the telluric ab- 
sorption features. Spectroscopic flat-fields were obtained 
from a halogen lamp installed at the Cassegrain port in 
front of the instrument's entrance window. The sky con- 
ditions were fair with the seeing between 0'.'45 and 0'.'75 
at R band. 

3. DATA REDUCTION 

The data were reduced using IRAF^ and custom- 
written IDL codes. The two dimensional spectrograms 
were pair-subtracted and then averaged and the pixel 
responses were normalized by the dark-subtracted flat- 
field images. Noisy pixels were filtered out based on the 
pixel statistics, and interpolated prior to spectral extrac- 
tion. The one-dimensional spectra were obtained using 
the IRAF aperture extraction package. 

The removal of the telluric lines by dividing the data 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



by the spectra of the standard stars was handled by IDL 
codes that take into account possible observational flaws, 
including slight shifts in the wavelengths, differences in 
airmasses (and therefore depths of telluric absorption 
lines), slightly different spectral resolutions, fringing of 
the continua, and saw-toothed features due to the differ- 
ent readout channels. In cases where the flaws were too 
extreme to be easily co rrected, a t ransmission spectrum 
calculated by ATRAN ()Lordlll99^ was used as a substi- 
tute to remove the telluric features. The signal-to-noise 
ratios of the spectra were calculated from the noises on 
the continua near the Hg absorption lines. 

Wavelength calibration was achieved by matching the 
observed spectra to the model transmission spectrum. 
The accuracy of the calibration is normally better than 
one tenth of one pixel (Sv — 0.75 km s~^) when the 
signal-to-noise ratio of the observed spectrum is suffi- 
ciently high. The radial velocities of the absorption 
lines observed in the multiple runs were corrected for 
the orbital motion of the Earth, converted to ulsr us- 
ing the radial velocity package of IRAF, and added to- 
gether with appropriate weights according to the signal- 
to-noise ratios of the spectra. The signal-to-noise ratios 
after adding up all available data range 60 to 300. 

The seven most relevant lines for this study are shown 
in Fig. [TJ For the Pistol Star, NHS 25 and possibly more 
weakly for NHS 22, the (3(1,0) line at 3.9530 ^^m is over- 
lapped at low wavelength with an atomic emission line. 
This wavelength region is omitted for NHS 25 in Fig. 1. 
The atomic emission matches the Si II lines at 3.9477 - 
3.949 9 /xm compiled in the NIST table (jRalchenko et all 
[2007I) . 

4. RESULTS 
4.1. Overview 

Ah five of the absorption lines (i?(l,0), i?(l,l)", 
i?(l,l)', (3(1,0), and Q{1,1)) that originate from the low- 
est para and ortho rotational levels, (J, K) — (1,1) and 
(1,0), were detected in all sight lines. In each sight line 
the line profiles of these transitions are basically the same 
(Fig- [J) 7 showing a broad trough of absorption primarily 
at negative radial velocities, and several sharp absorp- 
tion minima, which we attribute to intervenin g cloud s in 
the foreground spiral arms, as did lOka et al.l (|2005f ) for 
GCS 3-2. 

The 3.5337 /im R{3, 3)' line from the (3,3) metastable 
level was detected in all targets. In nearly every case 
it occurs largely at negative radial velocities. Its pro- 



Goto et al. 



TABLE 2 
Summary of observations. 



Source 


Exposure 


Grat 


ng" 


Line Coverage 








Atmo. 


ph 


eric Standard 


UT Date 




[s] 


ECH 


XDP 












Name 






Spe. 




GC IRS 21 . . . 


1800 


11050 


2800 


fl(2,2)' 


, -R(l,l)" 


, fl(l,0) 


Q(i,i) 


Q(i,o) 


HR 7528 






B9.5 IV 


2004 Jul 28 




2400 


5050 


3500 


fl(3,3)' 


, -R(i,i)' 


Q(i,o) 






HR 6378 






A2 V 


2004 Sep 02 


GC IRS 3 . . . 


360 


8350 


6100 


Rihiy 


, fl(i,o) 


«(i,i)' 


Q(1,1),Q(1,0) 


HR 7924 


(q 


Cyg) 


A2 lae 


2001 Jun 16 




1440 


4400 


5200 


i?(3,3)' 


, -R(i,i)' 








HR 7924 


(a 


Cyg) 


A2 lae 


2001 Jun 16 




3120 


11050 


2800 


fl(2,2)' 


-R(i,i)" 


, R{i,o) 


Q(i,i) 


Q(i,o) 


HR 7001 


(a 


Lyr) 


AO V 


2004 Jul 07 




2400 


5050 


3500 


fl(3,3)' 


-R(i,i)' 


Q(i,o) 






HR 7001 


(a 


Lyr) 


AO V 


2004 Jul 08 




3000 


11050 


2800 


fl(2,2)' 


-R(i,i)" 


, fl(i,o) 


Q(i,i) 


Q(i,o) 


HR 7121 






B2.5 V 


2004 Sep 25 




1280 


11050 


2800 


R{2,2)' 


-R(i,i)" 


, fl(i,o) 


Q(i,i) 


Q(i,o) 


HR 7121 






B2.5 V 


2004 Sep 26 


GC IRS IW . . 


2560 


5050 


3500 


fl(3,3)' 


-R(i,i)' 


Q(i,o) 






HR 7001 


(a 


Lyr) 


AO V 


2004 Jul 27 




2560 


11050 


2800 


R{2,2)' 


-R(i,i)" 


, fl(i,o) 


Q(i,i) 


Q(i,o) 


HR 7528 






B9.5 IV 


2004 Jul 27 


NHS 21 ... 


3600 


11050 


2800 


fl(2,2)' 


, i?(i,i)" 


: R{i,o} 


Q(i,i) 


Q(i,o) 


HR 7001 


(a 


Lyg) 


AO V 


2004 Jul 08 




3600 


5050 


3500 


fl(3,3)' 


, -R(i,i)' 


Q(1,0) 






HR 6378 






A2 V 


2004 Jul 08 


NHS 22 . . . 


1200 


11050 


2800 


fl(2.2)' 


, -R(i,i)" 


: fl(l.O) 


Q(i,i) 


Q(i,o) 


HR 7121, 


ATRAN 


B2.5 V 


2004 Sep 04 




1200 


5050 


3500 


fl(3.3)' 


, i?(i,i)' 


Q(1,0) 






HR 7121, 


ATRAN 


B2.5 V 


2004 Sep 04 


NHS 42 . . . 


1920 


11050 


2800 


R{2,2)' 


R{i,ir 


, RihO) 


Q(i,i) 


Q(i,o) 


HR 7924 


(q 


Cyg) 


A2 lae 


2004 Jul 27 




1920 


5050 


3500 


R{3,3)' 


fl(i,i)' 


Q(i,o) 














2004 Jul 27 


NHS 25 . . . 


1200 


5050 


3500 


R{3,3)' 


-R(i,i)' 


Q(i,o) 






HR 7121 






B2.5 V 


2004 Sep 04 




1200 


11050 


2800 


R{2,2)' 


-R(i,i)" 


, fl(i,o) 


Q(i,i) 


Q(i,o) 


HR 7121 






B2.5 V 


2004 Sep 04 


GCS 3-2 .. . 


120 


8350 


6100 


Rihir 


, fl(i,0) 


R{1,1)' 


Q(i,i),Q(i,o) 


HR 7528 






B9.5 IV 


2001 Jun 16 




240 


4400 


5200 


fl(3,3)' 


, -R(i,i)' 


Q(1,0) 






HR 7528 






B9.5 IV 


2001 Jun 16 




1344 


11050 


2800 


fl(2,2)' 


, R(i,ir 


, fl(i.o) 


Q(i,i) 


Q(i,o) 


HR 7001 


(a 


Lyr) 


AO V 


2004 Jul 07 




1344 


5050 


3500 


fl(3,3)' 


, -R(i,i)' 


Q(i,o) 






HR 7001 


(a 


Lyr) 


AO V 


2004 Jul 07 



file is broad and rather smooth, with a width ap- 
proaching 150 km s~^, and with no sharp features. 
Hg lines of such width have never been seen in the 
Galactic disk. Moreover, searches for the i?(3,3)' 
line in the disk along sight lines where the H;^ lines 
from the lowest two levels have clearly been observed 
jGcballc & Oka 1996; McCall ct al. 1998; 'Gcball e et aP 
[1999; McCall et al. 2002; Indriolo et al. 2007 ) have been 
unsuccessful ( Oka et al]|2005D . For these reasons and fol- 
lowing [O^^eF^I] (|2005f) . we attribute the entire i?(3,3)' 
absorption to the gas in the CMZ. The R{3, 3)' line pro- 
file mimics the broad trough in the J = 1 H^ absorption 
line profiles. We use this resemblance to differentiate be- 
tween the absorptions in the CMZ and those in the spiral 
arms. 

The i?(2, 2)' line at 3.6205 /xm was detected in one 
sight line only, that toward GC IRS 3, and only at ve- 
locities near -1-50 km s""'^. Because this is the first time 
that this line has been detected and because absorption 
by it is only 1 - 2 %, we have carefully examined its 
reality. Two telluric CH4 lines line at 2762.48 cm~^ 
(3.61994 fim) and 2761.36 cm"! (3.62141 /.tm) with peak 
absorption strengths of 20 ~ 30 % overlap the -R(2,2)' line; 
incomplete cancellation of them could cause absorption 
features near the observed feature. However, the wave- 
length match is not precise and moreover, the width of 
the i?(2,2)' hne (40 - 60 km s"^) is twice that of each 
of the telluric lines, which are nearly unresolved. As 
a further check we obtained independent observations 
of GC IRS 3 from the United Kingdom Infrared Tele- 
scope (UKIRT) in September 2004 and July 2005, using 
the echelle in the facility spectrograph CGS4. Although 



somewhat different in detail, the UKIRT observations 
consistently show a 1 - 2 % absorption dip at the same 
wavelength as the Subaru data. Thus the detection of 
this line toward GC IRS 3 is genuine. All observed spec- 
tra of the -R(2,2)' hne in GC IRS 3 are shown in Fig. [U 
together with the i?(3,3)' and i?(l,l)' absorption lines. 

In contrast to the Galactic disk where the detectable 
presence of H;^ is limited to relatively few sight lines and 
their narrow angular proximity, every sight line observed 
toward the CMZ has abundant H^ . Although the targets 
presented here are within 30 pc of Sgr A* in projected 
distance, the lengths of all of the absorbing columns of 
Kg are very large (significantly l arger than this se para- 
tion) as discussed below (see also lOka et al.l (|2005f )) and 
the observed velocity dispersion is very large, with ab- 
sorption occurring from —150 km s"^ to +20 km s^^. 
This range approximately matches the entire span of the 
radial velocity of the gas in the CMZ on the front-side 
of the center (fSawada et al.ll2004D . These characteristics 
suggest that the gas containing the observed R^ is not 
located in well-separated entities, but is continuous with 
a large volume filling factor. The large column densities 
of Kg in all sight lines are most obvious in the bottom 
traces of Fig. [1] where two of the strong H^ spectral lines 
from the (1,1) and (1,0) levels overlap and yield a blended 
feature of large equivalent width. 



4.2. 



at 



equivalent widths and column densities 



Kg lines from the J == 1 levels are composites of 
absorption in th e CMZ and in foreground spiral arms. 
lOka et al] (|2005f ) separated the two components of the 
-R(l,l)' absorption profile in GCS 3-2 by using the 2.3 fim 
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Fig. 1. — Spectra of Hg toward eight sources in the Central Molecular Zone. The _R(1,1)" and /?(1,0) absorptions in the bottom trace 
are separated by 35 km s~^ and the profiles overlap. The wavelength interval near Q(l,l) is heavily interfered by the telluric absorption of 
N2O, causing the low signal-to-noise ratios. The gap in the spectrum of NHS 25 in the Q(1,0) line corresponds to a strong atomic emission 
line, probably of Si II. 



V = 2 - i?(l) line of CO which is composed of sharp 
absorption features at the radial velocities (—50 km s~^, 
—30 km s~^ and km s~^) of foreground spiral arms with 
little or no CMZ component. The sharp features in the 
H^ absorption profile, which mimic the CO spectrum, 
were subtracted from the i?(l,l)' line and the remainder 
was a scribed to H^ in the CMZ (see Fig. 4 of lOka et al\ 
I2OO5I) . 

We found that this method was not effectively ap- 



plicable for the spectra in Fig. [T] because of the lower 
resolution of the IRCS (by a factor of 3) compared to 
the Phoenix Spectrometer of the Gemini South Observa- 
tory, and the low signal-to-noise ratios of spectra toward 
targets fainter than GCS 3-2. We used an alternative 
method for separation, in which we assume for each star 
that the velocity profile of the trough in the J — 1 lines 
is a scaled version of the _R(3,3)' profile. The method is 
depicted in Fig. [3l It crudely assumes uniform physical 
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Fig. 2.— Spectra of the R(2,2)' line at 3.6205 fim toward 
GC IRS 3 obtained at Subaru and UKIRT. The i?(3,3)' and /?(!,!)' 
hnes observed at Subaru are also shown. 



conditions in the absorbing clouds along the line of sight. 
We used it to separate the CMZ and spiral arm compo- 
nents of the (5(1,0) line as well as the i?(l,l)' line. We 
ignore the detailed velocity profiles of the observed spec- 
tra, and list integrated absorptions of the i?(l,l)' line 
and their ratios in the CMZ to the total (CMZ plus spi- 
ral arms). This ratio is also the ratio of colum n densities. 
We te sted the method for GCS 3-2 for which lOka et all 
(j2005|) used the high resolution spectrum obtained by 
Phoenix at the Gemini South Observatory and reported 
the ratio (7V(1, 1)gc/^(1, l)Totai = 0.68) by a different 
method of analysis using CO spectrum; the agreement is 
satisfactory. The interference of the (3(1,0) line by the 
atomic emission for NHS 25 mentioned in Section 3 does 
not seriously hamper the measurement since the scaling 
of the i?(3,3)' line to the (3(1,0) line is not much affected. 
The ratio of the equivalent width of the CMZ component 
to the total ranges from 0.5 to 0.9 for the eight sight lines. 
The sight line toward GC IRS 3 is exceptional in that 
it passes through two different categories of gas. In ad- 
dition to the gas with mostly negative velocity common 
to all eight targets, it shows absorption at -|-50 km s~^. 
This absorption is unique in two other ways: (1) it is 
the only gas in which the _R(2,2)' line is observable (in- 
dicating that the gas producing this velocity feature is 
higher density than the gas seen at other velocities) and 



(2) it is the only gas for which the column density of 
the (3,3) metastable level is definitely higher than in the 

(1.1) ground level. Therefore we consider the -1-50 km s~^ 
absorption separately. Measured equivalent widths of all 
spectral lines are listed in Table [3l For those transitions 
starting from the (1,1) level, both total equivalent widths 
and equivalent widths for the part of the absorption orig- 
inating in the CMZ are provided. 

The H;^ column densities in the (1,0), (1,1), (3,3), and 

(2.2) levels were calculated from the equivalent widths 
of the 0(1,0) R (l,lY, E(3,3)' , and E(2,2)' lines, re- 
spectively, as in iGeballe et al.l ()1999t ). from the equa- 
tion W\ = (87r^/3/iz/)A^ievei|A*P, using the theoretical line 
strengths |^p provid ed by J. K. G. Watson (also calcu - 
lable from the list of lNeale. Miller, fc Tennvson I [l996l ) . 
The results are given in Table [H The column density cor- 
responding to the 50 km s~^ absorption toward GC IRS 3 
is listed separately at the bottom of the table. The total 
column densities A^(H3 )totai are simply the sum of col- 
umn densities of the (1,1), (1,0), (3,3 ), and (2,2) levels 
since modeling of the thermalization l|Oka fc Epdr2004 ) 
shows that they are the only levels significantly popu- 
lated at the temperatures and densities of the region. 

For broad spectral lines such as those treated in this pa- 
per, uncertainties in the equivalent widths result mainly 
from the uncertainties in defining the continuum level 
rather than the standard deviations of individual data 
point. The listed uncertainties were estimated by multi- 
plying the standard deviation of the continuum level by 
the width of the relevant velocity range. 

4.3. Temperature and Density 

The column densities given in Table [4] were used to 
determine the temperatures, T, and densiti es, n, of the 
gas based on the steady state calculation bv lOka fc Epd 
(2004). Their diagram, in which the values of the pop- 
ulation ratios n(3,3)/n(l,l) and n(3,3)/n(2,2) are plot- 
ted as functions of T and n (their Fig. 4), is inverted 
in our Fig. 0] so that the observed population ratios can 
be used directly to determine T and n. The population 
rat ios toward GCS 3-2 from the velocity resolved data 
bv lOka et al.l (|2005f ) are included in the figure for com- 
parison. Most of the clouds in the Galactic center fall 
in the narrow parameter region of the T — n diagram, 
200 K< T < 300 K and n < 50-200 cm'^, with the 
sole exception being the -1-50 km s~^ component toward 
GC IRS 3 for which a higher temperature of ~ 400 K 
and a density of 200 cm^'^ were found. 

The upper limits to the number densities in Ta- 
ble [4] are based on the upper limits to the i?(2,2)' line 
equivalent widths. Lower limits may be roughly esti- 
mated fro ni the total e x tincti on of Ay ~ 30 mag re- 
ported by ICotera et al.l (|2000f) toward t he Quintuplet 
and (C entral Clusters and the estimate bv lWhittet et al.l 
(|1997f l that the visual extinction of the diffuse inter- 
stellar medium toward the GC is ~ 20 mag, that is, 
E{B ~ V) '^ 6.5. Ascribing all of the latter to the 
CMZ and using t he standard conversion factor from 
E{B-V) to N{K) (Bohlin. Savage, fc Drake 1[l978h gives 
iV(Il2) -- 2 X 10^2 cm-2. Assuming L = 100 pc yields 
n(Il2) ^ 60 cm^'^. This shows that the gas density is not 
much lower than the upper limit. Indeed the very exis- 
tence of observable Hjj" also argues against much lower 
densities, in which conditions it would be difficult for the 
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Fig. 3. — Scaling the i?(3,3)' absorption (top trace in each panel) to match the absorption trough of the R(l,l)' absorption (middle 
trace in black). The scaled /?(3,3)' is overlaid in red with ij(l,l)', and the difference spectrum is shown in the bottom trace. Vertical lines 
indicate the range of velocities where the profiles were scaled. For GC IRS 3, the scaling is done separately in the two velocity ranges, 
-200-20 km s'^ (red) and 20-80 km s"^ (blue). 



II2, which is needed to produce HJ, to exist. 

5. DISCUSSION 

5.1. The Ionization Rate 

An estimate of the ionization rate, C, in the Galactic 
center is the most important result of these Kg observa- 
tions. To avoid confusion, we use C(H) and C(H2) for the 
ionization rates of H and II2, respectively, in the Galactic 
disk where ionization is mainly due to cosmic-rays. We 
use C for the rate in the CMZ where the effect of X-rays 
and FUV radiation may be significant. C(H) and C(H2) 
are r elated approximately as C(H2) = 2^(11) (|Dalgarnol 
[20061) . 

Early studies of the cosmic ray ionization rate 
were made in order to understand the physical 
state of th e interstellar mediu m a nd its heating 
mechanisms (iHavakawa. Nishimura. fc Takava nagi 1961; 
Spitzer fc Tomaskol Il968t iField Goldsmith, fc Habinel 



1969( ). ISpitzer fc Tomaskol (|1968[ ) reported a minimum 
value of Q = 6.8 x 10~^^ s~^ from an extrapolation 
of the cosmic ray spectrum observed at Earth, and an 
upper limit of 1.2 x 10~^^ s~^ calculated from an as- 
sumed supernova frequency. Their Q is the total ion- 
ization rate for atomic hydrogen including the effect of 
secondary electrons, written here as C(H ). The rnost re - 
cent measurement of C(H), reported bv IWebben ()1998l ) 
using data from the Voyager and Pioneer spacecraft at 
distances up to 60 AU from the Sun, is C(H) > (3- 
4) X 10^^^ s~^, which should replace the lower limit 



used bv ISpitzer fc Tomaskol ()1968D . 

With the advent of the Copernicus satellite and its 
far-ultraviolet spectrometer, which allowed FUV stud- 
ies of ch emical abundances in t he diffuse ISM (for a 
review; ISpitzer fc Jenkinj Il975f ). column densities of 
molecules and radicals such as HD and OH could be 
measured. The reaction sequences leading to the cre- 
ation of these two molecules start from H+ and are thus 
useful for estimating C(H). Values of C(H) determined in 
this wav (lO'Donnell fc Watsonlll974l: iBlack fc Dalgarnol 



Il977t iHartguist. Dovle. fc DalgarnolTl978D are mostly a 
few times 10~^^ s~^, comparable to the lower limit given 
above. The co mprehensive model of the dif fuse interstel- 
lar medium bv lvan Dishoeck fc Blackl (|1986l ) suggested a 
somewhat higher value of C(H) « 7 x 10~ s'-^. How- 
ever they used an Hj^ dissociative recombination rate 
constant one thousand times lower than the current value 
which was popular in some circles at the time. They 
noted that if the thousand times higher constant were to 
be used, the value of C(H) would need to be increased by 
a factor of ^^ 10. 

The formation and destruction mechanisms for Hg are 
the simplest among all astrophysical probes of the ion- 
ization rate, and this makes Hg , when detectable, a 
powerful tool for measuring the ionization rate. While 
the R^ column densities observed in dense clouds, 
(0.4 - 2.3) X 10^^^ cm-2^ are consistent with C(H2) on 
the order of 3 x lO^^^ s'^ (Geball e fc Oka! fl996l : 
iMcCallet all 119990 . the surprisingly high A^(H+) ob- 
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TABLE 3 
Equivalent widths of absorption lines. 
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served in diffuse clouds, comparable to those in dense 
clouds in spite of the i r 10 times smaller vi s ual extinction 
(iMcCall et all fl998t iGeballe et all fl999l : iMcCall et al] 
|2002| ). require an order of magnitude higher C(H2)- Cru- 
cial to this conclusion are the laboratory measurements 
of the di ssociativ e recombinati on rate consta r it fcp. in 
Eg. ( 1) (lAmanol Il988: Larss on et al. I 119931: iLarssonl 
[2000t IMcCall e t al. 200 ^.120041). 

Tog ether, iMcCall et all (|2002f ) and llndriolo et all 
(|2007f ) have carried out a survey of Rf in the diffuse 
interstellar medium toward 29 nearby stars. Hg was de- 
tected in 14 of these sight lines with column densities 
of (0.6 - 6.5) X 10^"* cm~^. For these diffuse clouds the 
product (L is in the range (0.5 - 4.4) x 10^ cm s~^. Sep- 
arating C and L is not straightforward. Based on various 
methods of estimating L, which gave values of 2.2 - 31 pc 
for different clouds, the cosmic ray ionization rates were 
found to be in the range C(H2) = (1.2 - 7.4) x 10"^^ s'^ 
with an average of 5 xlO^^^ s~^ . Note that th e primary 
ionization rate of H, (p, used bv llndriolo et al. (2007) is 
relate d to C(H2) by 2.3Cp = ^fH^) (jGlassgold fc Langeil 
Il974[ ). This has established that the cosmic ray ioniza- 
tion rate in diffuse clouds is an order of magnitude higher 
than in dense clouds. This result might be understand- 
able in view of attenuati on of soft components of c os- 
mic rays in dense clouds ([Cravens fc Dalgarnolll978r i. if 
a large flux of low energy cosm ic rays is ubiquitous in 
the Galaxy (jMcCall et al.ll2003r ). or if higher energy cos- 
mic rays are trapped b y Alfven waves in diffuse clouds 
(jPadoan fc Scald l2005f ). The discrepancy between the 
high C(H2) and previous low C(H) values determined from 
HD and OH has been explained by Liszt (2003) as due 
to neutralization of H+ by charge exchange with small 
grains which had not been considered in the previous 
studies, in which the neutralization of H+ was assumed 
to be due to much slower radiative recombination with 
electrons. 

As seen in Table [3l the total column density of H;^ 
in the CMZ toward the eight observed stars ranges over 
(1.8 - 6.1) X 10^^ cm^^, more than an order of magni- 
tude higher than the column densities in diffuse clouds in 
the Galactic disk, suggesting even higher values of ( and 
L. Since we are interested in setting the minimum value 
of both C and L, we use the maxi mum value of / = 1 
and the minimum value oi Rx = 3 (Sodroski et al.|[l995l : 



(Ci)r 



= 7.4 X 



10-"cm3s-i7V(H+)totai. (2) 



lArimoto. Sofue. fc Tsujimoto1ll996[ ) in Eg. (1) and ob- 
tain the numerical expression 



To obtain Eg. (2), we have used fee ~ 7.7 x 10' 
calculated for T = 250 K from t he temperature de- 
pendence reported bv lMcCall et al.l (J200 4i) assuming the 
same electron temperature as the gas temperature, which 
may be justified in view of the fast thermalization of 
electrons with molecules and atoms. We used the C 
to H ratio in t he so lar vicinity of 1.6 x 10"'* given 
by ISofia et ahl (|2004D . which is close to the value of 
1.4 X 10~^, determined directly from infrared spec- 
tra of CO and H2 toward an embedded young star 
NGC 202 4 IRS 2 (iLacv et al.l [1994). and 1.3 x lO"'', 
found by lUsuda et al.l ( 20081 ) with the same technigue 
but toward NGC 7538 IRS 1 in the outer Galaxy. The 
observed iV(H|^)totai yield lower limits for (L, {(L)jnin 
= (1.3 - 4.5) X 10^ cm s~^. Although lower limits, the 
values of {CL)nun are 1000 tin ies higher than valu es of (L 
determined for dense clouds (McCall et al.|[l99 y) and 10 
time s higher than values f or diffuse clouds in the Galactic 
disk ([Indriolo et al.|[2007l ). We cannot cleanly separate C. 
and L based on current understanding of the CMZ, but 
discuss two extreme cases in the following sections. 

5.2. The Case of Large L in the CMZ 

lOka et all ([2005f) favored a large L and hence a high 
volume filling factor for the newly found warm and dif- 
fuse molecular gas. This was not only because the i?(3,3)' 
metastable lines with Large velocity widths and similar 
velocity profiles had been detected along all eight lines 
of sight toward sources scattered over 30 pc in projected 
distance from Sgr A*, but also because the deduced T 
and n of the gas were similar in all cases. The latter 
suggests that the absorbing diffuse material is contin- 
uous rather than contained in several separate entities. 
Further support for this interpretation is that the range 
of H;^ absorption velocities in each sight line spans the 
entire range velocities observed at radio and millimeter 
wavelengths with large apertures and interpreted as com- 
ing from the front side of the CMZ. Extending spec- 
troscopy of Hg over more widely spaced sight lines in 
the CMZ will constitute a more definitive test of the fill- 
ing factor. A long path length may also be favored over 
high ionization, in that there are limits to the value of 
C for efficient production of Hg as discussed in the next 
subsection. 

L is unlikely to be significantly larger than the radius of 
the CMZ, which for the region of interest is ~ 130 pc, as 
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estimated using Fig. lib of lSawada et al.l (|2004l ). Setting 
L = 100 pc, close to this maximum permissible value, 
the aforementioned range of values of (Ci)min yields the 
lower limits Cmin — (0.4 - 1.5) x 10^^^ s^^ for the various 
sight lines. For lower values of / and higher values of Rx , 
Cmm will be higher. 

The large extent of the warm and diffuse molecular gas 
postulated above conflicts with the previous concept of 
the gas in the CMZ, re presented, for example, in Fig. 9 of 
iLazio fc CordesI ()1998[ ) where the ultra-hot X-ray emit- 
ting plasma fills the CMZ. The warm and diffuse gas 
revealed in this paper cannot coexist with the ultra-hot 
plasma because the plasma electrons will be immediately 
cooled and molecules will be destroyed. Thus in the 
above model the volume filling factor of the ultra-hot 
gas would be reduced by a large factor. 

The large extent of the ultra-hot plasma gas 
in the CMZ was proposed because of observa- 
tions of extensive X-ray emission near the GC 
(e.g., [ Koyama et al.' 1989; Yamauchi & Kovama 1993; 
iKovama et al. 1996). Sunvaev. Markevitch fc Pavlinskyi 



gas with nearly 100% s urface filling factor toward the 
CMZ (see iGoldreich fc~Sridhar 2006, for a different in- 



Rodriguez-Fernandez. fc Martin-Pintadd 



(|1993f ) and iMarkevitch. Sunvaev. fc Pavlinskv I (|1993l ) 
observed similar X-ray emission, but interpreted it dif- 
ferently in view of the difficulty of confining such a 
hot gas by the gravitational potential of the Galac- 
tic center, and the extraordinarily large energy in- 
put required to maintain the high temperature. With 
the advent of the Chandra X-ray satellite, a great 
many X-ray point sources have been resolved in the 
GC. These have been interpr eted as explaining all 
(I Wang. Gotthelf. fc Lang 1 12002| ). or a significant frac- 
tion ( Muno et al.l 120031 ) . of the observed intense X-rays 
as being due to stellar sources such as cataclysmic vari- 
ables and young stellar objects. Additional observa- 
tional reports also argue against the ultra-hot plasma 
based on the close correlation between the stellar distri- 
bution an d the observed intensities of the X-rays i n the 
GC (e.g., [Warwick. Sakano. fc Decourchelle I I2006D and 
in the Galactic ridge (e.g.. lRevnivtsev et al.ll2006f ). How- 
ever, interpreta tion of th e X-ray emission remains con- 
troversial (e.g.. | Muno et al. 2004; Morris fc Navakshi^ 
[20061 IKovama et al.ll2007D . 

Yet another category of high temperature gas observed 
in the CMZ is the hot ga s observed in reco r ubina tion lines 
and other phenomena. iLazio fc CordesI (|1998f ) studied 
the A^-dependent hyper-strong radio- wave scattering and 
free- free radio emissions and absorptions and interpreted 
the data as evidence for Te ~ 10^ K and Ue ^ 10 cm"'^ 



terpre tation) . 

|200_5!) have observed fine structure 

tion lines toward a sample of 18 



and 



recombina- 
sources and re- 
ported even hig h er ele ctron densities of 30 - 100 cm~^. 
ILazio fc CordesI (|1998[ ) inferred that these regions were 
at interfaces between molecular clouds and the X-ray 
plasma and proposed (see their Fig. 9) that the high 
density molecular clouds with a volume filling factor 
of perhaps /v ~ 0.1 are surrounded by the scattering 
electron gas, and the rest of the space is filled with 
X-ray emitti ng plasma. Ou r observations and analy- 
sis (see also lOka et al.l I2005D c ontradict such a picture 
(jBoldvrev fc Yusef-Zadehll2006D . 

5.3. The Case of High ^ in the CMZ 

It is usually assumed that the density of high en- 
ergy {E > 100 MeV) cosmic rays is approximately uni- 
form over the Galaxy because of their long mean free 
paths. The uniform distribution on the sky of dif- 
fuse 7 ray emis s ion a ppears to lend support to this 
(e. g. iDigel et al.l l200lD. Recent detailed calculations 
bv Biisc hing et al.l ()2005f ) typically give a 20 % vari- 
ation depending on location and time. However, the 
situation is very different for low-energy cosmic rays, 
which are the most important for ionization of the gas; 
their energy densities can vary greatly dep ending on 
proximity to local cosrnic sources dKu lsrud fc Cesarskvl 
ll971l : ISpitzer fc Jenkinjll975HCesarskvlll975l ) because of 
their short mean free paths. For example, a proton of 

2 MeV, which effectively ionizes hydrogen, travels only 

3 pc befor e losing its energy in a me dium with uh = 
100 cm~3 (jCravens fc Dalgarnolll978[ ). The ionization 
rate by low-energy cosmic rays could be proportional to 
the number density of th e accelerating sources (SNRs; 
e.g. IKovama et al. I Il995f l and therefore have strong lo- 
cal variations. If the energy density of cosmic-rays fol- 
lows the surface density of OB stars, the cosmic ray 
ionization rate would be an order of magnitude higher 
in the Galact ic cent er th an in the solar neighborhood 
( Wolfire et all 120031). lYusef- Zadch. Law, fc Wardle I 
( 2002t l and lYusef-Zadeh et al.l (|2007. ) called for an even 
higher local enhancement of low-energy cosmic rays in 
the Galactic center (C = 2 x lO"^'* s'^ to 5 x lO'^^s-M) 
to account for the iron fluorescence line at 6.4 keV. Such 
values are much larger than those derived in the pre- 
vious section using H;^. X-rays and far ultraviolet ra- 
diation from stars may also increase C. A factor of 10 
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Fig. 4. — Temperature and densit y p l otted in filled circles as functions of population ratios n(3,3)/n(l,l) and n(3,3)/n(2,2) obtained 
by inverting the diagram of Oka & Egg ||2QQ3). Temperature is indicated by thick white line s and density by thin lines. The estimated 
uncertainties in n(3,3)/n(l,l) arc shown in error bars. Velocity resolved data of GCS 3-2 from [Oka et al.l II2005I ) is shown in filled squares 
for comparison. 



or higher increase of C due to X-ray emission has even 
been reported in t he circumsteUar disk of a low mass star 
(jDotv et al.ll2004l ). The enhancements could be large for 
the CMZ, where there are abundant intense emitters of 
X-rays and EUV radiation. 

For the above high values of C,, L would be quite short 
if Eq. (2) is applied. However, there are two arguments 
against adopting such high values of C- First, as argued 
earlier, the observed high velocity dispersion of Hg indi- 
cates that the gas is quite extensive rather than confined 
in many small clumps. This indicates that there are non- 
local sources of ionization covering an extensive region. 

Second, the results of Eqs. (1) and (2) are based on a 
simple linear formalism, which is a good approximation 
for low C- For much higher ( this should be replaced 
by higher order non-linear formalism in which A^(H3 ) no 
longer increases proportionally to C, but tends to satu- 
rate. One s ource of such nonlinearity, revealed by iLisztJ 
(|2006l[2007[) . is the dissociative recombination of Hj, H^ 
-|- e~ — > H -|- H, which competes with the Hj^ produc- 



tion reaction Hj 



Hs 



H^ 



H. The former reaction, 



in addition to the dissociative recombination of Hg , re- 
duces the H;^ concentration twice. It also increases n(H) 
and decreases /, hence the non-linea rity. In an e xample 
given in the lower figure of Fig. 3 of lLisztl (|2006[ ). a ten 
times increase of ( from 10~^^ s~^ to 10~^^ s~^ leads 
to only a fourfold increase of Hj^ (the saturation is more 
severe for the upper figure). The importance of satura- 
tion depends sensitively on the value of /, the fractional 



abundance of H2. Another source of non-linear behav- 
ior is the production of H through dissociation of H2 by 
cosmic rays. X-rays, and EUV radiation. The rates of 
these processes are comparable to the rate of ionization 
(Glassgold & Langcr 1973, 1974). Since the reproduc- 
tion of H2 on dust grains is slow, H accumulates in the gas 
and the charge exchange reaction, HJ -|- H — > H2 + H+, 
which has a high rate constant of 6.4 x 10"^*^ cm^ s~^, 
will compete with the Hg producing reaction introduc- 
ing additional non-linearity. The increase of n(H) and 
resulting reduction of / affect the gas in several other 
ways, all of which make the increase of ^ ineffective in 
increasing H^ . While more detailed treatments of these 
nonlinearities have yet to be worked out, it seems un- 
likely to us that C can be larger than '--^ 3 x 10~^^ s~^ 
over extensive portions of the CMZ. Small volumes with 
higher C may well exist in the CMZ, but they would not 
necessarily produce a large fraction of the observed H|^ . 

5.4. The +50 km s^^ absorption toward GC IRS 3 

The sight line toward GC IRS 3 is special in that at 
+50 km s^^ LSR it has a detectable population of Hj]" 
in the unstable (2,2) level in addition to the large popu- 
lations in the (1,1) and (3,3) levels. Also this is the only 
case for which the population in the (3,3) metastable 
level is definitely higher than in the (1,1) ground level 
(see A^(3,3)/7V(1,1) in Table 4). The observable pop- 
ulation in the (2,2) level indicates a high density and 
the larger population in the (3,3) level than in the (1,1) 
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level indicates high temperature environment. Our anal- 
ysis yields n = 175 - 300 cm'^ and T = 350 - 500 K. 
As seen from Fig. 4, this environment is clearly distinct 
from those elsewhere. We speculate below on two pos- 
sible locations for this gas: (1) the Circumnuclear Disk 
(CND) and, (2) the "50 km s^^ cloud" behind the GC. 

5.4.1. Circumnuclear disk 

GC IRS 3 is only 5" (0.2 pc) from Sgr A*, which is 
well within the projected span of the Central Cluster 
(the maximum transverse separation of bri ght s tars in 
the cluster is ~0.65 pc; iViehmann et al. 20051 ). The 
central 3 pc of the GC contains hot, high density gas 
which shows intricate mini-s piral structures (jEkers et al.l 
119831 iLo fc Claussml |1983[) termed the Western Arc, 
the Northern Arm, the Eastern Arm, and the Bar 
(jGiisten et al.l Il987l ). The extension of the Western 
Arc has been associated with a nearly completely cir- 
cular, Keplerian circumnuclear disk by analysis of HCN 
radio emission (iGiisten et al.l [19871 iGenzel fc Towned 
119871 iJackson et allll993t [Christ opher etld] l2005( l and 
by other methods; recent direct obser vations of the disk 
via it s far infrared dust emission by iLatvakoski. et aLl 
()1999( l are particularly noteworthy. There also 

have been attempts to interpret other three mini- 
spirals as CNDs which are not n ecessarily elliptical 
(iLiszt. Burton, van der Hulst 1 119851; [Q uinn fc Sussman] 
Il985t iLacv et al.l 119911 : IJackson et al.lll993t lLis^l2003L 
e. g.). GC IRS 3 is located in the cavity between the 
Northern Arm and the Bar, a region of low gas den- 
sity. H. S. Liszt (private communication) has suggested 
that the H I absorption at - 1-50 k m s~^ in Fig. 5 of 
iLiszt. Burton, van der Hulstl ()1985f ) may be due to the 
same cloud producing the Hj^ absorption in the sightline 
toward GC IRS 3. Not only does the velocity match, but 
also the high velocity gradient of 100 km s~^ arcmin"^ 
found bv ILiszt. Burton, van der Hulst I ()1985f ) is consis- 
tent with the width of the H^ 50 km s~^ absorption 
feature. 

However, there are two possible problems with this in- 
terpretation. First, it calls for an uncomfortably high 
cosmic ionization rate in order to account for the ob- 
served H^ column density in the short path length within 
the CND. Assuming a path length on the order of 0.5 pc, 
the observed (Ci)min of 7.5 x 10"* cm s^^ in Table[4]gives 
C ~ 2 X 10 ~^^ s~^, an extremel y high value. Such a value 
is possible (jYusef-Zadeh et aLll2007) in view of the high 
density of stars and radiation field in the re gion, but in 
such conditions the higher order chemistry (jLisztl l2006l . 
|2007() discussed in Section 5.3. surely plays a role. Al- 
though the detection of Hj^ in the unstable (2,2) level 
indicates higher density compared with other sight lines, 
the population of the (3,3) level is still much higher and 
the distribution is highly non-thermal. The density is 
definitely much lower than that reported for the CND 
(lO'^ - 10*^ cm-3). This could be due to the fact that 
GC IRS 3 is in the northern cavity and is only behind 
the edge of the CND rather than its densest regions. 

The second problem is the absence of the i?(2,2)' line 
in the spectrum of GC IRS IW. It is remarkable that 
the i?(2,2)' line is detected toward GC IRS 3, but not 
toward GC IRS IW which is jus t 8.^4 apart, or 0.31 pc 
for a common distance of 7.6 kpc (jEisenhauer et al1l2005l : 



iNishivama et al.ll2006r). If both s tars belong to the Cen - 
tral Cluster (jBecklin et~aI]ll978HViehmann et"al]l2005| ). 
it is not likely that they are much more than 0.5 pc apart. 
Although the source might well be located in the fore- 
ground, GC IRS IW overlaps with the high density part 
of the Northern Arm, while GC IRS 3 is located in the 
cavity. It is difhcult To understand why only GC IRS 3 
shows i?(2,2)' absorption, if the absorption occurs within 
the CND. 

5.4.2. The"50 km s-^ cloud" 

The observed velocity of -1-50 km s"^ toward GC IRS3 
immediately suggests that the absorption might arise 
in the well known "50 km s~^ cloud", one of the 
giant molecular clouds (> IO^Mq) located near the 
Galactic center. The "50 km s~^ cloud" is some- 
times referred to as M— 0.02— .07, o r was called the 
"40 km s-i cloud" b y lOortJ (119771). and has been 
discussed in detail by iBrown fc Lisztl ()1984D . Many 
atomic and molecular studies have established that 
it is located behind the Galactic nucleus, in be- 
tween the non-thermal radio source Sgr A East (SNR) 
and the H II region Sgr A West, which enshrouds 
the central star cluster fWhitcoak. R ogstad. fc Lockhart 
igTi"; 'Glisten fc Hcnkell983; Pauls, J ohnston, fc WilsoU 
1996) . Yusef-Za deh et al. (1996) argued that the several 



masers they discovered in the southeast of Sgr A* are ex- 
cited in the shock-heated interface between Sgr A East 
and M— 0.02— 0.07, where the expan ding supernova rem - 
nant hits the giant molecular cloud. iTsuboi et al.l ()2006l ) 
also found clues for an interaction between Sgr A East 
and the "50 km s^^ cloud" in the high velocity wings 
of CS J = 1 — 0, although it may not be the case that 
the "50 km s~^ cloud" is the sole cloud that falls in the 
radial velocity -1- 50 km s~^ in the line of sight toward the 
G alactic center ([Zvlka. Mezger. fc Lesch 1992). 

iGeballe. Baas, fc Wadd (|1989D found that GC IRS 3 
and IRS 7 show deep absorption lines of the CO funda- 
mental band at 4.7 fim at the radial velocity -1-50 km s~^, 
while GC IRS 1 and IRS 2 do not. They attributed the 
+50 km s"""^ absorption to the CND but also mentioned 
the possibility that GC IRS 3 is not a member of the 
central cluster, but behind the -1-50 km s~^ cloud and 
thus not in the Central Cluster. The nature of GC IRS 3 
is controversial (|Pott et al.ll2008^ ; it is among the bright- 
est sources in the Central Cluster region only less lum- 
nous than GC IRS 1 and GC IRS 10, and much redder 
than them in the mid-infrared even though it is (appar- 
ently) in the cavity. Because of its location and fea - 
tureless infrared continuum (e.g. iKrabbe et al.|[l995[) . 
GC IRS 3 has widely been believed to be a hot mas- 
sive star heavily obscured by its own dust (similar to 
the bright Quintuplet sources) , in addition to foreground 
dust. However, new interferometric observation that re- 
sol ye the circu r tistell ar dust shell have been interpreted 
by iPott et all ()2008l ) as evidence that the central ob- 
ject is a late-type carbon star, the photospheric spectral 
features of which are completely obscured. This inter- 
pretation is consistent with the apparent lack of ioniza- 
tion of the circumstellar shell, which one might expect 
if the object lies within the Central Cl uster, and indeed 
the la ck of ionized gas near GC IRS 3 (jRoche fc AitkenI 
|1985() . Although somewhat unlikely due to its proxim- 
ity on the plane of the sky to the Central Cluster, the 
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placement of GC IRS 3 well behind the cluster and the 
"+50 km s~^ cloud" is otherwise a straightforward as- 
signment and avoids the difficulties associated with the 
CND hypothesis mentioned in the previous section. 

5.5. Infrared pumping 
Durin g the 2005 Royal Society Discussion Meeting 



on H^ ([O^] [20063) , J- H. Black proposed an infrared 
pumping mechanism which can produce a non-thermal 
rotational level distribution of H^ if the infrared flux is 
sufficiently high. The mechanism is illustrated in Fig. [5l 
Infrared radiation pumps Hg from the fully populated 
lowest ortho-level, (1,0), to the (2,0) level in the first 
excited vibrational state via the i?(l,0) transition. The 
molecule decays to the (3,0) level in the ground state 
by spontaneous emission with a lifetime of 25.3 ms. 
The molecule further decays to the (3,3) metastable 
level through the forbidden rotational transition (3,0) 
— > (3,3) e mitting a far-infrared photon at 49.62 /im in 
3.75 h ours (|Pan L Okalll986HNeale. Miller, fc Tennyson I 
Il996f ). Thus, HJ molecules are pumped from the (1,0) 
level to the metastable (3,3) level. 

Since the spontaneous emissions are all very fast, 
the rate determining process for Black's mechanism is 
the rate of infrared pumping of the (1,0) level of H;^. 
We rewrite the excitation rate given in Eq. (15) of 
Ivan Dishoeck fc Blackl (|1982[ ) as 



Wr^ia) — A 



?7r= 1.51 X 10 



-7^ 



where (f),^ is the photon flux of the ambient radiation at 
the wavelength of the transition in the units of photon 
cm~^ s~^ Hz~^, and the values of A = 3.6685 /xm, the 
Einstein coefficient A = 98.5 s~^, and the dilution fac- 
tor of 0.286 are used. The dilution factor takes into ac- 
count the spontaneous emission (2,0) — > (1,0) back to the 
ground level (1,1) with A = 35.5 s~^, resulting in a reduc- 
tion of the pumping efhciency. In order for pumping to 
be significant, the value of W_r(i.o) has to be comparable 
to the collision rate, 



R = nav = nk^ 



2 X 10~'^s~\ 



where we used a number density of n = 100 cm ^ and 
the temperature independent Langevin rate constant, 
2 X 10"^ cm^ s-i for the H;j" + H2 colfision. 

The energy density of the interstellar radiation field in 
the solar n eighbor hood at 4 /im ha s been estimated by 
iPorter. Moskalenko. fc Strong! (flxjg ) to be Aua = 0.055 
/im eV cm~'^ fjxaT^ which translates to (/ij^ = 6 xlO^^ 
photons cm~ ^ s~^ Hz~^. This is consistent with t he value 
reported by iMathis. Mezger. fc Panagial (119831 ) within 
a factor of 1.5. The corresponding pumping rate is 
Wji^ifi) = 1 X 10^^-^ s~^, about 2x10^ times lower than 
the collision rate. iPorter. Moskalenko. fc Strong! (J2006! ) 
have estimated an average energy density about 20 times 
higher in the GG. Launhardt, Zylka, & Mezger (2002) 
have reported an energy density in the central 120 pc 
(their Fig. 10) which is another factor of 10 higher. The 
corresponding pumping rate is still 100 times lower than 
the collision rate and the infrared pumping is perhaps 
not very effective for majority of the gas reported in this 
paper. 

In the central parsec of t he Galaxy, however , the stellar 
density is ~ 10^ M© pc'^ (jSchodel et al.!!2007f) . including 
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Fig. 5. — Schematic of the infrared pumping mechanism pro- 
posed by J. H. Black. H3 in the lowest ortho-level (1,0) is excited 
to the (2,0) level in the first excited vibrational state via the ij(l,0) 
transition (red line) by ambient photons. The excited H^ sponta- 
neously decays to the (3,0) ground level through the P(3,0) tran- 
sition (orange line) and then to the metastable (3,3) level (shown 
in blue). See text for more details. Other transitions used in this 
paper are also shown. 



a large number of lumi nous red giants and a t least 80 
hot and massive stars (jPaumard et al.! I2006D . In this 
environment infrared pumping likely would overwhelm 
the collisional relaxation of any H;^ that is present. The 
effect, however, must be local to the very central region 
of the GMZ and probably does not seriously alter the 
population of the (3,3) level beyond a few tens of parsecs 
from the center, i.e., in the bulk of the GMZ. However, it 
may well be the best explanation for the strong i?(3,3)' 
absorption line toward GC IRS 3 at 50 km s^^ discussed 
in the previous subsection. 

6. SUMMARY POINTS 

(1) Seven infrared spectral fines of II|^, (5(1,0), (5(1,1), 
i?(l,l)', i?(l,0), R{l,lY, i?(2,2)', and i?(3,3)', have been 
observed toward eight bright and hot stars distributed 
from the Quintuplet Cluster ~ 30 pc east of the Galactic 
Center, to the Central Cluster close to Sgr A*, using the 
IRGS spectrograph of the Subaru Telescope. 

(2) All sight lines showed strong i?(3,3)' absorptions 
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with high velocity dispersions, indicating a high popu- 
lation of the (3,3) metastable rotational level. Those 
absorptions are ascribed to H;^ in the CMZ and indicate 
a large surface filling factor. 

(3) The five absorption lines starting from the low- 
est levels, (1,1) and (1,0), are composed of sharp spec- 
tral lines due to H;^ in the foreground spiral arms and a 
broad trough due to H^ in the CMZ. These two absorp- 
tion components have been separated by assuming that 
the profile of the trough is proportional to the i?(3,3)' 
absorption profile. 

(4) The observed high H^ column densities indicate 
lower limits of (L, (C-Zj)min = (1-3 - 4.5) x 10^ cm s~\ 
which are 1000 and 10 times higher than for dense and 
diffuse clouds in the Galactic disk, respectively, indicat- 
ing high values of C and L in the CMZ. 

(5) Various arguments suggest C is on the order of 
10"^^ s"i and L on the order of 50 pc in the CMZ, 
indicating a high ionization rate in the region and a high 
volume filling factor of the gas. 

(6) Analysis based on the thermalization model calcu- 
lation (jOka fc E pd 2004) has given T = 200 - 300 K and 
n < 50 - 200 cm^^ (Fig.H]). A lower limit on the order of 
50 cm~'^ is set for n from the observed visual extinction. 
Ay, and the existence of detectable Hj. 

(7) Taken as a whole the observations reveal a previ- 
ously unidentified gaseous environment in the CMZ hav- 
ing a high volume filling factor. This is radically differ- 
ent than the previous picture of the gas in the CMZ, and 
suggests in particular that the filling factor of ultra-hot 
X-ray emitting plasma may be nowhere close to unity. 

8) The sight line toward GC IRS 3, close to Sgr A *, 
is unique in having a detectable population in the Hg 
(2,2) unstable level near -1-50 km s~^, indicating a higher 
density for the gas producing this absorption. We specu- 
late that the cloud responsible for this absorption feature 
is either a compact cloud near the circumnuclcar disk or 
the well known "50 km s^^ cloud" . For the former a high 



value of C = 2 X 10^^** s^^ is needed, while for the latter 
GC IRS 3 is located behind the "50 km s~^ cloud" and 
thus outside of the Central Cluster. 

(9) An infrared pumping mechanism to populate the 
(3,3) metastable level proposed by J. H. Black has been 
examined. The infrared radiation in the extended CMZ 
is not intense enough to significantly pump this level, and 
the high temperatures detected there are likely thermal. 
However in the CND pumping probably dominates and 
temperatures determined for gas in that region based on 
the relative populations in this level and the J = 1 levels 
would not be kinetic temperatures. 
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